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Abstract 

The storage and retrieval of a single-photon polarization q-bit by 
means of STIRAP through the atoms with degenerate levels is studied 
theoretically for arbitrary polarization of the driving laser field and ar- 
bitrary values of the angular momenta of resonant atomic levels. The 
dependence of the probability of long-term photon storage on the polar- 
ization of the driving field and on the initial atomic state is examined. 



1 Introduction 

The search for an efficient optical quantum memory and its implementation for 
quantum networks is an active area (see, e.g., the reviews [T} [HI 131 H] ) - The role 
of carrier of information - flying q-bit - in such devices is played by a single 
photon and the most natural way to encode the q-bit is provided by its two 
polarization degrees of freedom (see [5j[6l[7l[8] for some recent proposals). While 
in the majority of the proposed schemes the photon quantum state is mapped 
to the quantum state of some atomic ensemble [SJ [7J |S] , in the experiment [5] 
the photon polarization q-bit was recorded in the superposition of magnetic 
substates of the long-lived degenerate level of a single rubidium atom. Such 
single-atom quantum memory offers an advantage of processing of the stored 
quantum information. Since at the storage stage the state of a single atom, in 
which the photon q-bit is encoded, is well determined, it may be altered in a 
controlled way, e.g., by some laser pulse, so that the state of the retrieved photon 
will differ from that of the recorded one in a desirable way implementing thus 
the performance of quantum gates. The setup employed in [5] for the photon 
storage and retrieval was based on the same technique of vacuum stimulated 
Raman scattering involving adiabatic passage (STIRAP), which was previously 
successfully employed for deterministic single photon emission [9J 1101 111] . In 
these experiments the three-level A - type atom was trapped inside the high- 
finesse cavity, one branch of the atomic A - type transitions was coupled to the 
quantized cavity field, while the other one was coupled to the driving coherent 
laser field. The STIRAP with the three-level atom with non-degenerate levels 
is well described in the reviews and textbooks (see, e.g., [HI H2]). In case of 
non-degenerate levels there exists the only dark state - the superposition of the 
ground state a and some metastable state b - uncoupled to the excited level c 
(Figure 1). In course of Raman scattering the atom is adiabatically transferred 
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Figure 1: The level diagram 



from the initial state b to the target state a. The STIRAP with degenerate 
atomic levels and with classical coherent resonant fields was studied in [T3] , 
were it was shown that unlike the non-degenerate case the population from 
the initial state b is not always totally transferred to the target level a. The 
objective of the present paper is to analyze the process of storage and retrieval 
of a single-photon polarization q-bit by means of STIRAP through the atoms 
with degenerate levels for arbitrary values of the angular momenta of resonant 
atomic levels and to study the dependence of the probability of long-term photon 
storage on the polarization of the driving field and on the initial atomic state. 

In section 2 the interaction model is described and the instantaneous eigen- 
vectors of the interaction operator, which determine the evolution operator 
in the adiabatic approximation, are constructed. In case of degenerate levels 
there appear the new types of these eigenvectors, non-existing in case of non- 
degenerate levels, like the dark states, which atomic part belongs to only one of 
the lower levels a or b, and the bright states, which couple the excited level c 
with only one of the lower levels a or b. In section 3 the formula for the proba- 
bility of long-term storage of the photon polarization q-bit is obtained and the 
conditions for the photon storage with unit probability independent on its po- 
larization are outlined. The transitions with the angular momenta Jj, = J c = 1, 
J a = 2, corresponding to the transitions between the hyperfine structure com- 
ponents of the electronic levels 5 2 Si/2 and 5 2 Pi/2 of the 87 Kb atom, which were 
employed in the experiments [5], are analyzed. 

2 Basic equations 

We consider the three-level A - type atom in a high-finesse cavity, one branch - 
the transition J a — > J c with the frequency uj c q is coupled to the driving coherent 
laser field, while the other - the transition Jb — > J c with the frequency ujq is 
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coupled to the quantized cavity field (Figure 1), J a , Jb and J c being the values 
of the angular momenta of the levels. The electric field strength of the coherent 
laser field may be put down as follows: 

E c = e c (t)l c e- iUat + c.c, 

while the quantized field of the cavity in the interaction representation is de- 
scribed by the operator: 

E = e(i)(a 1 li+a 2 l 2 )e- iwt + /i.c., 

where the carrier frequencies lo c and w are in resonance with the frequencies w c o 
and ujq . e c (t) and l c are the slowly varying amplitude and the unit polarization 
vector of the driving field, e(t) is the slowly varying amplitude of the photon 
field, li and 1 2 are the two unit orthogonal vectors of the two polarization 
modes of the photon field, d\ and ci 2 are the photon annihilation operators 
for this modes. The temporal dependence e{t) of the cavity field amplitude 
appears due to some tailored alteration of cavity parameters. The equation for 
the slowly-varying density matrix p of the system, which consists of a single 
three-level atom and two-mode cavity field, in the rotating-wave approximation 
and in case of Raman resonance uj c q — co c = loq — lj = A is as follows: 



dt P 2 



V(t),P , (1) 



V(t) = -2AP c + G(i) + G t (t), ( 2 ) 
G(t) =n a (t)g a + Q b {t)g b , (3) 

9b = 9bia[ +56 2 4- ( 4 ) 

Here P c is the projector on the subspace of the atomic excited level J c , fl a (t) = 
2\d a \e c (t)/H and Qb(t) — 2\db\e(t)/h are the reduced Rabi frequencies for the 
coherent laser field and for the cavity field, d a = d(J a J c ) and c4 = d(JbJ c ) 
are the reduced matrix elements of the electric dipole moment operator for the 
transitions J a — > J c and Jb — > J c , while 

9a = g a l*, 9u = g b l*, i = 1,2, 

g a and g b are the dimensionless electric dipole moment operators for the tran- 
sitions J a — > J c and Jb — > J c . The matrix elements of the circular components 
9aq (a — a,b; q = 0, ±1) of these vector operators are expressed through Wigner 
3J-symbols [TT)] : 

(g aq )rn a ,m c = (-l) J ^ ( J « 1 Jc V (5) 

The solution of the equation (flj is expressed through the evolution operator 
S(t): 

pit) = S{t)p{Q)S + {t). (6) 
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In the adiabatic approximation the evolution operator S(t) is defined by the 
instantaneous eigenvectors \vk(t) > and eigenvalues Afc(i) of the interaction 
operator V(t) in a following way: 

S(t) =5^exp{i^ fc (t)}|u fc (t) >< w fc (0)|, (7) 
k 

Mt) = \ J Mt'W . (8) 

Since only single photon storage and retrieval are discussed in the present paper 
and no decay or dephasing processes are taken into account it is sufficient to 
limit the system space to the subspace with the basis vectors \ J a m a > |0,0 >, 
\Jbm h > 1 1, >, \ Jbm b > |0, 1 > and \ J c m c > |0,0 >, where \ J a m a >, \ Jbm b > 
and | J c m c > denote the atomic Zeeman states, while \ni,n 2 > (rii,2 = 0, 1) are 
the field number states with n\ photons in the first polarization mode and n 2 
- in the second. This subspace with the dimension N = 2{J a + 2J b + J c + 2) 
constitutes the invariant subspace of the interaction operator V(t), so that its 
matrix represents itself a square hermitian N x N matrix. The states \a > 
|0, > and |c > |0,0 >, which atomic part belongs to the level a or c, may be 
represented as columns with 2 J a + 1 or 2 J c + 1 elements correspondingly, while 
the states |6i > 1 1, > +|6 2 > |0, 1 >, which atomic part belongs to the level b, 
may be represented as columns with 2(2Jb + I) elements: 

|6i > |1,0> +|6 2 > |0,1 > = 

Then the operator g a will be represented by the (2J a + 1) x (2J C + I) matrix, 
while the operator g b will be represented by the 2(2Jfc + 1) x (2J C + 1) matrix 

were each block gu (i = 1, 2) represents itself a (2 J b + 1) x (2 J c + 1) matrix. 

In order to find out the instantaneous eigenvectors \vk(t) > and eigenvalues 
Xk{t) of the interaction operator V(t) let us start with the operator 

&(t)G(t) = nl(t)glg a + Ql(t)gtg b , 

which acts at the subspace of the upper atomic level c (note that g\gb = 0, 
g\g a = 0). The eigenvectors of this operator form an orthonormal set, while 
its eigenvalues are real and non-negative. The eigenvectors \D^(t) > with zero 
eigenvalues: 

&(t)G(t)\D c n (t)>=0, n = I,...,iV c d , (9) 

uncoupled to the lower levels a and 6, are at the same time the eigenvectors of 
the interaction operator V(t) with the eigenvalues \ cn — — 2A. The eigenvectors 
\C n (t) > with positive eigenvalues c^(t) > 0: 

&(t)G(t)\C n (t) >= c 2 n (t)\C n (t) >, n=l,...,JV C) (10) 
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are coupled by the electric dipole transitions to the states 

\F n (t) >= -±-G{t)\C n {t) > . (11) 

Cn{t) 

These states \F n {t) > also form an orthonormal subset at the subspace of the 
lower atomic levels a and b, as it follows from (fT0|) - (ITT1) . 

With the introduction of the states \C n {t) > and \F n {t) > all the eigenvec- 
tors of the interaction operator V(t) with non-zero eigenvalues may be easily 



obtained. This operator has 2N C eigenvectors 

\VW(t) >=sm6 n (t)\F n (t)+ cos 6 n (t)\C n (t) >, (12) 

\Vy\t) >= cos e n (t)\F n (t) >-nmO n (t)\C n (t) >, (13) 

tan2^(t)=-^, (14) 

with the eigenvalues 

A( ± )(<) = -A± v^+4(f). (15) 



All the states \D^(t) > and (V^^i) > constitute the orthonormal set of = 
2(2 J c + 1) — N d eigenvectors of the operator V(t) with non-zero eigenvalues. 

The other N d = N — eigenvectors \D n {t) > of the interaction operator 
V(t) obtain zero eigenvalues. These states - dark states - with the atomic part 
belonging to the subspace of the lower atomic levels a and b remain uncoupled 
to the upper level c: 

&(t)\D n {t) >=0, n = l,...,N d . (16) 

As it follows from (|TT1) and (|16p . all the dark states \D n (t) > are orthogonal to 
all the bright states \F n (t) >. 

Among all the dark states let us distinguish first the dark states \D% > and 
|D„ >, which atomic part belongs to the lower level a or b only. These states 
are time independent and may be obtained as the eigenvectors of operators ffoffl 
or gi,gl with zero eigenvalues: 

9a9l\D a n >=0, k = l,...,N d , (17) 

g b gl\D b n >=0, k = l,...,Nl (18) 

The other N d b eigenvectors |-D° b (£) > of the interaction operator V(t), which 
atomic part belongs to both atomic lower levels a and b, obtain zero eigenvalues 
and satisfy the equation (|16|) . while 

gt\D£{t) >^0, gl\D^{t) >^ 0. 

The temporal dependence of these dark states may be immediately obtained 
from the equation (|16|) : 

\D*(t) >= z n (t)[n a (t)\B d n > -n b (t)\A d n >], 
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where Z n (t) is the normalization factor, while \A d l > and \B„ > are temporally 
independent states, which atomic parts belong to the levels a and b correspond- 
ingly and which satisfy the equation 

gl\Ai>=g\\B d n >^Q. (19) 

Introducing the matrix 



D b = ^ w \B n ><B n 



b 2 

n n 

containing only the eigenvectors \B n > of matrix g b gl with non-zero eigenvalues 
6^, we may write the equation (fT9|) as follows: 

\S* >= D ba \Ai >, D ba = D b g b gl (20) 

Now we may define the ortho normal set of states \A~ > as the eigenvectors of 
the hermitian matrix b\ a D ba with non-zero eigenvalues: 

Dl a D ba \A d n >= a 2 dn \A d n >, a dn > 0. (21) 
Then, as it follows from (|2U|) and (|2T|) , the states 



\Bi >= —t) ba \A d n > 



(22) 

O-dn 

also constitute the orthonormal set, so that the orthonormal set of dark states 
\Dn b (t) > ma y be expressed through these states \A d > and \B d > by the 
formula: 

1 n (t) Tm^um ■ (23) 

All the eigenvectors of the interaction operator V(t), comprising the set of states 
\D c n {t) >, \V^ ±] {t) >, and the set of dark states \D^ >, \D b n > and \Df{t) > 
constitute the complete orthonormal set of states, which determines the evolu- 
tion operator ([7]) in the adiabatic approximation. 

3 Photon storage and retrieval 

Initially the atom is at the lower level b, its state being defined by the atomic 
density matrix /)q, while the state of the field is a pure single-photon state 

|/>=Ci|i,o>+6|o,i> ! |£i| 2 + I6I 2 = i, 

where the two complex numbers £i and £2 define the photon polarization (the 
polarization q-bit), so that the initial density matrix /5(0) of the atom-(- field 
system is: 

/ 3(0)=pM/x/|. 



G 



The classical coherent laser field is adiabatically switched off, while the interac- 
tion with the quantum field is adiabatically switched on in the time interval T\ , 
so that: 

0,(0) = n oll n„(Ti) = o, n 6 (o) = o, n 6 (ii) = n w . 

The atomic part of the initial bright states \F n (0) >, defined by the equation 
(fTTj) . with non-zero eigenvalues c^(0) ^ belongs to the level a only, so that the 
corresponding initial eigenvectors | (0) > of the interaction operator, defined 
by (p~2]) ~ (|14p . will not contribute to the evolution. However the atomic part 
of the initial bright states |-F ra (0) > with zero eigenvalues c^(0) = 0, obtained 
as the limit at t —> +0, belongs to the level b. Since c^(0) = ^ al c^ n , where 
(? an are the eigenvalues of the operator g\g a , such states will be coupled to the 
eigenvectors |C* > of the operator g\g a with zero eigenvalues: 

gtg a \C b >=Q, n = l,...,N b c . (24) 

So, the evolution operator |(7J| will be determined by the normalized states 
g b \C b n > and the dark states \D b n > and \Df(t) >, defined by (JISJ) and (23), 
which atomic part also belongs to the level b. The states > in the evolu- 
tion operator are responsible for the photon passage through the cavity without 
interaction, while the states <?b|C„ > are responsible for the process of photon 
absorption and the adiabatic transfer of the atom from the initial level b to the 
excited level c to some state, which is uncoupled to the lower level a. In both 
cases the final system state is unstable and for purposes of long-term photon 
storage the presence of such states in the initial system state should be avoided 
by the proper choice of polarization of the driving field and by the preparation 
of the atomic state. With such states excluded the evolution is determined only 
by the states \D? b (0) >, which are responsible for the long-term recording of 
the photon polarization q-bit onto the superposition of Zeeman substates of the 
metastable level a. The evolution operator at the instant of time t = T\ after 
the recording process is finished is as follows: 

^ = -Y^K >< B d n \, (25) 

a 

so that each state \B!~ > from the level b, defined by (f2"2"|) . is adiabatically 
transferred to the corresponding state \A% > from the level a, defined by (f2"Tj) . 
The total probability w of the photon storage at the level a generally depends 
on the photon polarization £i and £2: 

2 

w = tri^msl} = J2 "V ( 26 ) 

where 

<%• = E < b n \Po\ b n >, (27) 

n 

while Ib^ 1 >=< 1, 0\B d > and \bf? >=< 0, > are the pure atomic states at 
the level b - the projections of the field+atom states |B„ > on the two photon 
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polarization states 1 1 , > and |0,1 >. For purposes of q-bit storage the total 
probability must be unity and must not depend on photon polarization, so that 
the condition 

Wij = 5 l . ] (28) 

must be satisfied. 

To retrieve the photon after the storage time r the classical coherent laser 
field with the same polarization l c , as that of the field applied to store the 
photon, is adiabatically switched on, while the interaction with the quantum 
field is adiabatically switched off in the time interval T 2l so that: 

fia(Ti + r) = 0, n a (Ti + r + T 2 ) = Q a2 , 

n 6 (Ti + t) = n b2 , n 6 (Ti + r + r 2 ) = o. 

If the photon was stored with unit probability, then at the start of storage 
process the atom was prepared at some state of the level b, which contained only 
the parts \B% > of the "dark" states \D^ b >, while the the initially unpopulated 
level a was populated during the storage process by the adiabatic transfer of 
the parts \B% > of the "dark" states \D® b > to their counterparts \A* > at 
the level o, so that at the start of the retrieval stage only the "dark" states 
|D° h > will be present in the initial atomic state, provided that the polarization 
of the retrieving field is the same as that of the recording field, ensuring that 
the "dark" states \D® b > are the same for both stages - the storage and the 
retrieval. In this case only the dark states \D^(t) > will contribute to the 
evolution operator ([7]). Then, at the end of the retrieval process at the instant 
of time t 2 = T\ + r + T 2 the evolution operator 

S2 = -Y.\ Bd n>< Ad n\ ( 29 ) 

n 

is determined by the same states |A„ > and \B% > as the evolution operator Si 
(|23|) and the product of the evolution operators S 2 Si is reduced to the projector 
on the subspace of states \B% >: 

n 

and since only such states are present at the initial system state /5(0), then this 
state will be retrieved without distortion 

h = S^mslsl = p(0), (30) 

reproducing a single photon \f > with the recorded polarization and the atom 
at the initial state with the density matrix p\. 

Let us choose the two unit polarization vectors li and 1 2 of the quantum 
field in the plane XY as the two vectors with opposite circular components: 

hq = Sq,-ll hq = $q,l, 
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while the polarization vector l c of the classical field generally contains also Z- 
projections. 

In the experiment [5] the angular momenta of the atomic levels were J a = 2, 
Jh = J c = 1, while the driving field was linearly polarized along the axis Z: 
Icq = $q,o (7T-polarized field). In this case there are three states \A% > at the 
level a, defined by ([2"T]) : 



\A( >= 








0) 















> \At >= 


I 














W 








w 



which are adiabatically coupled to the three states \B* > at the level b, defined 
by p2]l: 

'0 s 

|1,0>, \B d 2 >=- I 1 1 |0,1>, 




\&i >= -7= 
3 V2 



The elements of matrix u 

IT 

Pq in a following way: 




1,0 > - 




0,1 > 



i,j <|27|) determining the photon storage probability 
are expressed through the elements p^ b m , of the initial atomic density matrix 



W 2 .2 



bb , i bb 
- n bfc 4- I„ bb 



W2,l =W 12 



■-n bb 
2 Pl,-l, 



so that the conditions f|28[) may be met only if the atom is initially in the pure 
state \Jb — l,m(, = >, as it was in the experiment. Then, as it follows from 
the expression (1251) for the evolution operator, the initial system pure state 

\J b = l,m b = 0> (£i|l,0 > +6|0, 1 >) 

is adiabatically transferred to the pure state 

(&| Ja = 2,m a = -1 > -&|J„ = 2,m a = 1 >)|0,0 > . 



The conditions (|28p may be also met, as calculations show, with the driving 
field linearly polarized along the axis X: 

l cq = - S q s) 
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and with the same initial atomic state |J& = 1,7715 = >, though in this case 
the photon polarization will be mapped to the superposition of states \ J a — 
2, m a — > and | J a = 2, m a — ±2 >: 

a°|m a = > +a + |m a = 2 > +a~|ra a = -2 >, 

gi+jg ± V3(6+6) 1 

a = a = ± . 

2V2 4 2 

For the lower values Jb — and Jb = 1/2 of the angular momenta of the initial 
atomic level b the conditions (|2"8|) for the photon storage with unit probability, 
independent on photon polarization, may be easily met. With the driving field 
linearly polarized along the axis Z: l cq = 6 q fi, these conditions will be satisfied, 
for example, for the transitions with Jb = 0, J a — J c = 1 and with Jb = 1/2, 
J a = Jc = 3/2 for any initial atomic state, in particular, for the equilibrium 
state: p h \ /2 _ 1/2 = P?/ 2 ,i/ 2 = V 2 , P-i/2,1/2 = pf/2,-1/2 = °- However for 
larger values of angular momentum Jb > 1 it becomes hard (if not impossible) 
to meet the conditions (|28l) . Indeed, the memory efficiency (or better to say 
inefficiency) is determined by the number N\ of the "dark" states of the type 
\D b n >, defined by the equation (fT8|) , and by the weight of these states in the 
initial atomic state. The number of states 2(2J{, + 1) with the atomic part at 
the level b is doubled due to the two possible polarizations of the cavity field 
and only 2 J c + 1 of them are coupled to the excited level c, while the other 
N b d = 2(2J b - J c + 1) - 1 states - the "dark" states of the type \D b n > - stay 
uncoupled, reducing the memory efficiency. The number of such states obtains 
its minimum N% — 2Jb — 1 on the transitions J5 — > J c = J5 + 1, so that at 
Jb > 1 the presence of these states becomes mandatory. The number N% of 
the states \D„ > grows with Jb and at Jb > 1 it becomes highly improbable 
that each Zeeman state 1 75,7715 > will not be contained in these "dark" states 
\D b n >. The general approach, developed in the present paper, enables one 
to analyze numerically the memory efficiency for any reasonable values of the 
angular momenta of resonant atomic levels. For example, let us consider the 
transitions Jb = 2 — > J c = 3 — > J a = 4 and the 7r-polarized driving field 
(Icq = Sg t o). Then the elements of the probability matrix Wij (|2T)) are expressed 
through the elements p b ^ b m , of the initial atomic density matrix pQ in a following 
way: 



1 66 1 1 66 1 6 55 55 5ft . . 

= =P-2,-2 + 2^-1 "I + jPofi + Pl,l + ^2,2> (31) 

— 1 n bb 4- 1 n bb 4- 6 n bb 4- n bb 4- n bb fW\ 

W2,2 - ljP2,2 + 2^1,1 + jPOM + P-l-l + P-2,-2, (32J 



^66 



while Wx,2 an d W2,i are expressed only through the non-diagonal elements p mb 
(nib 7^ fn'b) of the atomic density matrix. In addition the populations of Zeeman 
sublevels are non-negative p b ^ lm > and the total population of the level is unity 

A-2 + A-i + <o + Pm + P2.2 = I- (33) 
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For a hundred per cent memory efficiency w±.i — W2.2 — 1, then extracting (|3ip 
from (|33|) we obtain 

6 bb , 1 hb 1 1 6fc n 

jP-2-2 + 2^-1,-1 + 7^0,0 = 0, 

which is true only if p hb 2 ,_2 = = Pao = 0, an d extracting from (j3"3"j) 

we also obtain that p\ 2 = p b \i = 0, so that in this case the photon q-bit cannot 
be stored with unit probability. However, if the atom is initially prepared at 
the state | = 2, m b = >, then lUy = (6/7)5^, and the photon will be stored 
with the probability w = 6/7 independent of its polarization. 

4 Conclusions 

In the present paper we have obtained the general expression for the probabil- 
ity of long-term storage of a single-photon polarization q-bit at the magnetic 
substates of some long-lived degenerate atomic level by means of STIRAP. The 
dependence of this probability on the polarization of the driving field and on 
the initial atomic state was studied. It was shown that the conditions for the 
photon long-term storage with unit probability independent on the photon po- 
larization may be satisfied with low values of angular momentum of the initial 
atomic level Jb = 0, 1/2, 1. At larger values of this angular momentum J& > 1 it 
becomes hard (if not impossible) to store the arbitrarily polarized photon with 
unit probability. The general approach, developed in the present paper, enables 
one to analyze numerically the memory efficiency for any reasonable values of 
the angular momenta of resonant atomic levels. It was also shown that, since the 
conditions of the photon storage with unit probability are satisfied, the driving 
field with the same polarization as that used to store the photon will retrieve the 
initial atomic state and a single photon with the recorded polarization without 
distortion. 
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